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In this chapter we consider perturbations and stabihty of higher dimensional black holes 
focusing on the static background case. We first review a gauge-invariant formalism for linear 

■ perturbations in a fairly generic class of (m-l-n)-dimensional spacetimes with a warped prod- 
' ' uct metric, including black hole geometry. We classify perturbations of such a background 



introduce a set of manifestly gauge invariant variables. We then introduce harmonic tensors 
' _| ' ' and write down the equations of motion for the expansion coefficients of the gauge invariant 

"•pi , perturbation variables in terms of the harmonics. In particular, for the tensor-type pertur- 

I ■ bations a single master equation is obtained in the (m -I- n)-dimensional background, which 

is applicable for perturbation analysis of not only static black holes but also some class of 
^ ' rotating black holes as well as black-branes. For the vector and scalar type, we derive a set 

, , of decoupled master equations when the background is a (2 -I- n)-dimensional static black 

hole in the Einstein-Maxwell theory with a cosmological constant. As an application of the 
' master equations, we review the stability analysis of higher dimensional charged static black 

^ \ holes with a cosmological constant. We also briefly review the recent results of a generaliza- 

, tion of the perturbation formulae presented here and stability analysis to static black holes 

' in generic Lovelock theory. 
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§1. Introduction 

There is a large variety of black hole solutions in higher dimensions. The stability 
of such exact solutions is clearly an important issue, and analysing linear perturba- 
tions of the existing exact solutions would be the first step to take. If a stationary 
black hole solution is shown to be stable under perturbations, it implies that the 
solution describes a possible final state of dynamical evolution of a gravitating sys- 
tem. If, on the other hand, an instability is found, it then indicates the existence of 
a different branch of solutions, which the original solution may decay into, and one 
can then anticipate more variety of black hole solutions. 

Apart from the stability issue, perturbation analysis also tells us a lot about 
basic properties of black hole solutions. The spectra of quasinormal modcs^^ contain 
information about the geometric structure of the background metric, especially near 
the horizon. The study of stationary perturbations of a stationary black hole solu- 
tion provides a criterion for the uniqueness/non- uniqueness property of the solution. 
Considering such stationary perturbations of a known solution may also be useful 
when attempting to construct approximate solutions. 

In this article, we shall review the perturbation formulae for higher dimensional 
static black holes and the stability analysis, following, to a large extent, Refs. 2)- 
4). The linearised Einstein equations off of a black hole spacetime are in general 
still quite involved, having a number of perturbation variables intricately coupled, 
and one therefore needs to simplify them to a tractable form. In 4-dimensions, 
the perturbed Einstein equations for a stationary vacuum black hole solution can 
be reduced to a set of simple decoupled ordinary differential equations (ODEs) by 
exploiting some particular geometric feature of the background solution: the Regge- 
Wheeler-Zerilli equations^) '^^ for the Schwarzschild metric and Teukolsky equations^^ 
for the Kerr metric. For higher dimensional rotating black holes, we are still a long 
way from having such a complete formulation for perturbations, though considerable 
progress along this direction has recently been made for some special cases^^"^^) [see 
Chapter 7 and references therein for the rotating black hole case]. Fortunately, 
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for static black holes, e.g., Schwarzscliild-Tangherlini metric and its cousins, such a 
reduction is possible in arbitrary higher dimensions, and a set of decoupled master 
equations, which correspond to the Regge- Wheeler- Zerilli equations in 4-dimensions, 
are now available.^)' More precisely, consider a (2 + n)-dimensional static black 
hole in Einstein-Maxwell system with a cosmological constant as our background, 
in which the n-dimensional internal space is an Einstein space and describes the 
horizon cross-section geometry. Then perturbation variables are classified into three 
types according to their tensorial behavior on the internal space. For each type of 
perturbations, a basis of gauge-invariant variables is introduced, and the Einstein 
and Maxwell equations are written in terms of them. For each type of perturbations, 
the perturbed equations of motion are further reduced to a set of decoupled master 
equations for a single scalar variable on the 2-dimensional part of the background 
spacetimc. Furthermore, the master equations thus obtained in Refs. 2), 3) take 
the form of a second-order self-adjoint ODE with respect to the radial coordinate 
of the black hole background and are therefore immediately applied to the stability 
analysis.^)' 

In the next section we first establish our notation and convention by describing 
our background geometry. We next explain how to decompose tensor fields in our 

background spacetime and how to construct manifestly gauge-invariant perturbation 
variables. In section 3, we introduce harmonic tensors on n-dimensional Einstein 
space and provide several theorems concerning basic properties of the harmonic ten- 
sors. We also give explicit expressions for scalar, vector and tensor harmonics on 
S"" in terms of the homogeneous coordinates for S"" in Subsequently from 

section 4 to section 6, we focus on perturbations of static black holes in the Einstein- 
Maxwell system with cosmological constant and derive a set of master equations. In 
section 7, we describe how to derive the master equations for static black holes in 
most general Lovelock theory recently obtained. Then in section 8, using the mas- 
ter equations, we show that a large class of static black holes in the Einstein-Maxwell 
system with cosmological constant are stable with respect to linear gravitational as 
well as electromagnetic perturbations. We also briefly comment on the recent study 
of instability of Lovelock black holes. Section 9 is devoted to summary and discus- 
sion. 

§2. Notation and conventions 

In this section we describe our background spacetime and discuss how to classify 
tensor fields on the spacetime. We then discuss the problem of gauge freedom and in- 
troduce gauge-invariant variables. We generally follow the notation and conventions 
of the papers 2), 4), and 3), which throughout this chapter we refer to as Papers I, 
II, and III, respectively. 

2.1. Background Geometry 

We consider an (m -|- n)-dimensional spacetime whose manifold structure is lo- 
cally a warped product type, ^ = x J^", and accordingly we often need to 
distinguish between tensors living in these different manifolds, ,A^"^, and .J^". 
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For this reason, we do not employ the abstract index notation^^^ in this chapter, 
and instead, we use upper case latin indices in the range K, L, M, N, . . . to denote 
tensors on lower case latin indices in the range a,b, . . . ,h on ,yV"^, and lower case 
latin indices in the range i, j, . . . ,p on J^". Accordingly, we introduce coordinates 
= {jj°',z^) in terms of which our background metric is written 

gMNdx^dx^ = gab{y)dy''dy'' + r\y)j{z)ijdz'dz^ . (2-1) 

We assume that an m-dimensional metric gab{y) on .yT"* is Lorentzian and an n- 
dimensional internal metric 'y{z)ij on J^'" is Einstein i.e., 

Rij = {n- l)Kjij (2-2) 

for some constant K, with Rij being the Ricci tensor of 7jj. When is maximally 
symmetric, the constant K corresponds to the sectional curvature of J^", and in 
what follows we normalize K = 0, ±1. We assume that J^" be complete, as it 
describes the geometry of cross-sections of the event horizon. 

We denote by Vm, Da, and Di, the covariant derivative operators compatible 
with qmn, dabi and 7jj, respectively. Having these derivative operators, we can define 
the curvature tensors on and .J^", and find their relations in terms of the 

coordinate components as 

R°'bcd = "^R"'bcd , (2-3a) 

R\bj = :^9^3 , (2-3b) 

R'jkl = R'jkl - {Drf{5i^ji - Sjjjk) , (2-3c) 

where '^R"'bcd and R^jki are the curvature tensors of gab and jij, respectively. Prom 
this and Eq. (2-2), we obtain 

(2.4a) 
(2.4b) 



Rab— 


^"'Rgab-n 


Rai~ 


0, 


Rij = 


[-- + (" 


R = 


Or 

""R - 2n— 



9i3 , (2-4c) 
K - {Drf 



(2.4d) 



where □ = D"-Da and "^R is the scalar curvature of gab- Note that the Ricci tensor 
takes the same form as in the case in which Jf^" is maximally symmetric. 

The geometric structure of our background spacetime requires that the back- 
ground stress-energy tensor Tmn should take the form 

T„, = 0, T'j = P5'j, (2-5) 

where P is a scalar field on 
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Then, the background Einstein equations, including cosmological constant, A, 
are written as 

DaDbr ( n{n - I) K - {Drf Ur\ 

_ y„.-.j,,.-.j..-y^,j ^ _ 1) - ^ _ ^ ^ (2-7) 



"^ab - - ' ' - n— gab = -Agab + i^'T^b , (2-6) 

r \ 2 r^^ r J 

Imp (n - l}{n - 2) K - (D r)^ , ^^Dr 
2 2 r2 



where denotes the gravitational constant. Note that although our main concern 
is about static black holes for which m = 2 with the metric form given by Eq. (2-14) 
below, our metric ansatz above also allows us to consider various different geometries, 
such as black-string/branes when m > 3 and Myers-Perry black holes with a single 
rotation when m = 4. 

2.2. Static black holes in (2 + n)- dimensions 

Now, having static black hole geometry in mind, let us set m = 2 and consider 
an electromagnetic field ^mn as a source for the background gravitational field. 
The field strength ^mn may be given by 

^ = ^Eoeabdy'' A dy'' + ^^ijdz' A dz^ . (2-8) 

Then, from V[m=^jvl] = 0, we obtain 

Eo = Eo{y), ^ij = ^ij{z), a[fe^y]=0, (2-9) 
and from Vn-^^^ = 0, 

= Viv=^"^ = ^e'^'Dbir^Eo) , (2-lOa) 

= V N-^'^ = Dj^'^ ■ (2- 10b) 



These equations imply that the electric field Eq takes the Coulomb form, 

E^ = ^^ (2-11) 

and # = \^ij{z)dz^ A dz^ is a harmonic form on Jif^. Although, in general, there 
may exist such a harmonic form that produces an energy-momentum tensor consis- 
tent with the structure of the Ricci tensors in Eq. (2-4), in the following we consider 
only the case J^ij = 0. 

With this assumption, the energy-momentum tensor for the electromagnetic 
field, 

= ^ML^J- - ^gMN^LK^'-'' , (2-12) 

is written 

T(-^)a^ = -PS\, T^^^^'j = PS'j; P=^El = ^. (2-13) 
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Then, solving the background Einstein equations, we have, when Vr / 0, the black 
hole type solution 

ds^ = -f{r)dt^ + 7(0 ^ '"''^'"^ ' ^^'^^^ 

where 

/W = i^-Ar^-^ + ^, (2-15) 

and 

A:=^^, Q^:=-^. (2-16) 

The spacctimc described by this metric can contain a regular black hole for some 
restricted ranges of the parameters, M, Q, A, K. [For the allowed parameter regions 
for regular black holes, see Appendix A of Paper III.] 

Note that we can also consider, as our background, solutions with Vr = 0, 
for which ^ becomes a cartesian product of a 2-dimensional maximally symmetric 
spacetime =yK^ and the Einstein space Jf"'. Such solutions include the Nariai solu- 
tion,-'^^)'-'^^) as well as the Bertotti-Robinson solution in 4-dimensions. The solutions, 
and their parameter ranges are given in Paper III and also in Ref. 17). As in Pa- 
per III, we can study perturbations of such Nariai-type solutions in a similar manner 
as in the black hole background case, but in this chapter, we are not going to deal 
with this case. 

2.3. Decomposition of vectors and symmetric tensors on compact manifolds 

When considering perturbations, we will in general have to deal with two major 
issues. One is concerning the ambiguity in a choice of gauge, due to the invariance 
of the Einstein equations under an infinitesimal gauge transformation. This issue 
will be discussed in the next subsection. The other one is that even when linearised, 
the Einstein equations are still a set of intricately coupled equations for a number 
of perturbation variables of the metric and matter fields, and are in general very 
difficult to solve. We therefore need to reduce the linearised Einstein equations to a 
simple tractable form. For this purpose we first classify perturbation variables into 
three different types according to their tensorial behavior on in such a way that 
the linearised Einstein equations get decoupled and can be dealt with separately for 
each type of perturbations. We then introduce harmonic tensors on so that 
each type of the perturbed Einstein equations reduces to a set of equations on the 
m-dimensional spacetime {^'^,gab)- This procedure enables us to obtain a set of 
significantly simplified equations for master scalar variables as we will see in later 
sections. 

Let us consider how to classify perturbation variables. We first note the following 
two decomposition theorems: 

(i) Suppose (^",7ij) be a compact Riemannian manifold. Any dual vector field on 
J^'^ can be uniquely decomposed as 



Vi = Vi + DiS 



(2-17) 
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where D^Vi = 0. This is essentially the well-known Hodge decomposition theorem, 
and we refer to Vi and S, respectively, as the vector- and scalar-type components of 
the dual vector Vi. 

(ii) Suppose {J(f"^,jij) be a compact Riemannian Einstein space, Rij = cjij for some 
constant c. Any second rank symmetric tensor field tij can be uniquely decomposed 
as 

tij = t'if + 2D^it^^\) + tL^ij + LijtT , (2-18) 
Lij := DiDj - -7yA (2-19) 

where IJ^tJJ^ = 0, = q, D¥^\ = 0, and ti = t™m/n. We refer to t^f , \ and 
{tT,tL), respectively, as the tensor-, vector- and scalar-type components of tij. Note 

(2) 

that tensor component t^j exists only when n > 3. 

A similar decomposition theorem in which is considered to be maximally 
symmctric-has been proved by Kodama and Sasaki. A general proof of the above 
theorems (i) and (ii) is found in Ref. 19), in which the compactness of J^'"' is used 
in an essential way to show that any symmetric elliptic operator, such as &D(^iVj^ 
on Vi, is essentially sclf-adjoint and the spectrum of any self-adjoint operator is 
discrete. The symmetric, elliptic operators appeared in the proof will be essentially 
self-adjoint on many non-compact manifolds of interests, even though their spectra 
will be continuous, and analogous decomposition theorems should also hold for many 
non-compact manifolds. In the following when considering the case in which is 
non-compact, we simply assume that analogous decomposition results hold. 

Now let us consider metric perturbations Hmn = Sqmn on our background 
spacetime {^,gMN)- We may project Hmn relative to the Einstein manifold 
as 

hMNdx^dx^ = habdy^dy^ + 2haidy''dz' + hijdz'dz^ . (2-20) 

The component hab is purely scalar with respect to transformations on As for 
the components hai and hij , applying the above decomposition theorems (i) (ii) , we 
further decompose them into their scalar, vector and tensor parts with respect to 
JT" as 

hai = Diha + /li? , (2-21) 
hij = hPij + 2^(i4^^j) + hLlij + LijhP , (2-22) 

where 

D^h^ij = hP\ = 0, (2-23) 
L»«/ii^^=0, &hPi = 0. (2-24) 

Thus, the tensor part of /imat is h^ij, the vector part of hMN consists of {h^^^ ,hj:\), 
and the scalar part of hj^N consists of {hab, ha,hL, hj!^). 
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Similarly, we can decompose perturbations of the energy-momentum tensor, 
STmn, into the tensor part 5Tj?)j, the vector part {6T^l\ STj!^), and the scalar 
part {STab, STa, STl, STt), where 

5Tai = DidTa + 5T^^ , (2-25) 
STij = 6T^^. + 2D^,5Ti^^j) + Sn-fij + LijdTr , (2-26) 

with b^5T!^\j = 0, Sr'^'^h = 0, D^5T^^^ = = D'6T^p,. 

The linearised Einstein equations are decomposed into these three types, and 
thus we can deal with perturbations of each type separately. Before writing down 
the equations for each type of perturbations, we discuss the gauge issue. 

2.4. Gauge-invariant formulation 

As is well-known, the Einstein equations are invariant under gauge transforma- 
tion generated by an (infinitesimal) vector field . Accordingly, a perturbation 
variable and its gauge transformed one, such as Sqmn and SguN — &^gMN^ both 
must describe the same physical situation, hence giving rise to an ambiguity in the 
representation of perturbation variables. In order to remove this gauge ambiguity 
and extract the physical degrees of freedom, we may proceed either by imposing ap- 
propriate conditions that completely fix the gauge freedom, or by constructing man- 
ifestly gauge-invariant variables and writing down the relevant equations in terms of 
them. The two approaches are equivalent. To see this, let us have a look at gauge 
transformation law of perturbation variables. 

In terms of the coordinates (y"',z^), the gauge transformation of the metric 
perturbation Hmn = ^9mn is written as 

hab ^ hab - Da^b " ^feCa , (2-27) 
hai ^ hai - r^Da (^^^ - A^a , (2-28) 

hij ^ hij - 2D^i^j) - 2r{D''r)iali3 ■ (2-29) 

Similarly, the gauge transformation of perturbation of the energy-momentum tensor 
is given by 

5Tab ^ 5T^b - eDcTab - TacDbC - TbcDae , (2-30) 
5Tai ^ 5Tai - Tabbii"" - r^PDa{^i/r^) , (2-31) 
STij ^ 6Tij - eDa{r''P)^ij - P{biCj + Dj^i) . (2-32) 

We can of course classify the gauge transformations above into the three ten- 
sorial types, and discuss each type separately. As is clear from the decomposition 
theorem (i), the generator has only the vector- and scalar-type component; it is 
decomposed as 

^a = Ta, ^i = Vi + biS, (2-33) 

where ^"^^DiVj = 0. Any tensor-type perturbation variable is by itself gauge-invariant: 
we have {hj:\j, 5T^\j). For the vector type perturbations, the gauge transformation 
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law of the metric perturbation is given as follows: 

hS^hS-r^D.{^, (2.34) 

h?i ^ - ■ (2-35) 

Inspecting the transformation laws, we find the following combination is gauge- 
invariant: 

i^W-fcll'-^oi^) (2-36) 

Similarly, inspecting the gauge-transformation law of the matter perturbation, we 
find two gauge-invariant variables: 

rl^ ■■= St!^ - , (2-37) 

T^^ := 2^(,<54'),) - 2Pb^ih^^\) . (2-38) 

Any vector-type gauge invariant variable can be expressed as a linear combination 

i^ai^ ' "^aP' '^ij^^ their derivatives. 
The above equations, (2-36), (2-37), (2-38), defining the gauge invariant variables 

may be viewed in a way that {h^^^ , ST^l^ , 6T^\) are expressed in terms of the gauge- 
invariants {Fai^ jT^KTip), and h^j. Since, under gauge-transformation, h^^ behaves 
just like — ^j, it can be chosen to take any value, and therefore one may view that 
h^j alone is responsible for the gauge ambiguity. This, in turn, implies that the 

specification of h^j in terms of {^ai^ ■, t-j^ ) , corresponds to fixing the gauge 
freedom of the vector type perturbation. 

For the scalar-type metric perturbation, the gauge transformation law is explic- 
itly given as follows: 

hab ^ hab - 2D(J'^^ , (2-39) 

ha^ha-Ta- r^Da , (2-40) 

hL^hL- 2r{D\)Ta - ^AS , (2-41) 
hT^hT-2S. (2-42) 

Now let us define Xm = {Xa, Xi = DiXi) by 



Xa := -ha + "^Da fe) , Xl := . (2-43) 



Then, noting that Xm gauge-transforms as Xm Xm + Cm, i-e., 

{Xa,XL)^{Xa + Ta,XL + S), (2-44) 
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we can immediately find gauge-invariant combinations: 

F^'^Kb = hab + 2L'(aX6) , (2-45) 

= hL + 2r{D\)Xa + -AXl . (2-46) 
n 

As for the matter perturbation, we find gauge-invariant combinations: 

:= STab + X'^DcTab + TacDbX'' + TbcDaX^ , (2-47) 

:= b,5T, + T„bA^' + r^PDa j > (2-48) 

:= (5Tl -PhL + r'^X'^DaP , (2-49) 
77(0) i^- _ ^.^.^j.^ + 2PLi^Xi . (2-50) 

Any scalar-type gauge invariant variable can be expressed as a linear combination 
of the gauge-invariant variables F^, i:(o)„6, T^, TT^jj) and their 

derivatives. 

One may view the above defining equations, (2-43), (2-45), (2-46), (2-47), (2-48), 
(2-49), (2-50), in such a way that all the scalar-type perturbations {hab, ha, fiL, hr) 
and {6Tab,5Ta,STL,dTT) are expressed in terms oi F^^\b, F^^\ U^^Kb, S^^Ki, ^(°\ 
i7(°)ij-which are gauge invariant, and -which can be taken as completely arbi- 
trary as the generator is. Then, one may say that specifying Xm in terms of 
the above set of the gauge-invariant variables (or assigning Xm some specific value) 
corresponds to fixing the gauge freedom. For example, the specification, Xm = 0, 
corresponds to the longitudinal gauge often used in the cosmological context, when 
m = 1, and to the Regge- Wheeler gauge, when m = 2. 

The perturbed Einstein equations are gauge-invariant and therefore can be 
written in terms of the gauge-invariant variables introduced above. Note that at 
this point all variables are functions on being dependent upon the (m -|- n)- 
coordinates, x'^ . In the next section, we introduce tensor harmonics on so that 
by expanding the perturbation variables in terms of the tensor harmonics, we can 
separate variables and reduce the relevant equations in ^ to equations in 

As a specific example of a source for gravitational field, let us consider the 
Maxwell field in the m = 2 case. Perturbation of the field strength 6^mn satisfies 
the perturbed Maxwell equations: 

y[M8^NL]=Q, (5 (ViV=^^^) = , (2-51) 

where here and in the following the external current is treated as a first-order 
quantity. The first equation implies that 5^mn is expressed in terms of the per- 
turbation of the vector potential, 5^/m, as S^mn = 2V^m^-^n]- Therefore 6^mn 
does not contain any tensor-type perturbation. The second equation gives two set 
of equations, 

i^A(r"(5^)"^) + bii6^r + ^oe"' l^lDbih\ - h%) - Wi'b^ = r ,(2-52a) 
1 



({6^Y + E^e'^Xb)] + DjS^i' = Ji , (2-52b) 
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where note that (6^)^^ = g^^g^^d^LK- 

The contribution of electromagnetic field perturbations to the energy-momentum 
tensor are given by 

^Tir^ = ^ {e'^S^cd + Eoh%) gab - ^E^ohab, (2-53a) 
^2-(em)a. ^ - Eoe'''' 6 .^i , (2-53b) 
STi<>^)i. ^ _^ (e'^H^cd + EQhi) d'j. (2-53c) 

For vector- type perturbations, we note that S^mn is gauge-invariant as well as 
invariant under a coordinate gauge transformation — >■ + . We find 

^^(em)(i) ^ _E,ea,DH4''^ , Sr^;"^^^'^ = , (2-54) 

where D'^S.^/^^^ = 0. Wc find from the first of Eqs. (2-52b) that a vector perturbation 

of the Maxwell field can be expressed in terms of the gauge-invariant variable Sj^^^^ 
on ^y^^ as, 

S^ab = , S^ai = Da6j^^'^ , 6^ij = DiSj^j'^ - D^l^^ . (2-55) 
For scalar-type perturbations, S^mn gauge-transforms as: 

S^ab ^ S.^ab - D^iE^e^ab , ^^aj ^ 5^a, " EoeabDjC' ■ (2-56) 

Note that 5,^ij = for a scalar perturbation. So, we can immediately find gauge- 
invariant combinations, 'S'a'^\ given by 

e„6<f = 6^ab + eabD,{EoX') , (2-57) 
eabDiS^'^^'' = + EoEabDiX'' . (2-58) 

The gauge-invariant variables introduced above are then written as 

^(^-)(°^a6 = -f {F('\b - F('K9ab) - Eo^^'^gab , (2-59) 
^(cm)(0)^. ^ .^o^.^o) ^ (2-60) 

^(em)(0) ^ ^2 ^^^^(0) _ ^F^O)c^^ , (2-61) 

Jjiemm,, ^ Q _ (2-62) 

§3. Harmonic tensors on the Einstein space 

When we write down perturbation equations and solve them, it is often more 
convenient to expand perturbation variables into Fourier-type harmonic components 
in terms of harmonic tensors appropriate for each tensorial type on the internal 
space J^". In this section, we summarize the definitions and basic properties of 
such harmonic tensors relevant to the descriptions in the present paper. 
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3.1. Scalar harmonics 

3.1.1. Definition and properties 

Covariant derivative D along J(^^ appears only in the form of the Laplace- 
Beltrami operator in the perturbation equations for scalar-type variables because 
'jij is the only non-trivial symmetric tensor on the Einstein space Hence, the 

scalar-type perturbation equations reduce to a set of PDEs on by expanding 
scalar-type perturbation variables in terms of scalar harmonic functions on J^T"' that 
satisfy 

(a + A;^)s = 0. (3-1) 

Here, when J^T'^ is non-compact, we assume that —A is extended to a non-negative 
self-adjoint operator in the L^-space of functions on J^f"^. Hence, k'^ > 0. Such an 
extension is unique if is complete^'') and is given by the Friedrichs self-adjoint 
extension of the symmetric and non-negative operator — A on C^{J(^^). 

If is closed, the spectrum of A is completely discrete, each eigenvalue has 
a finite multiplicity, and the lowest eigenvalue is k^=0, whose eigenfunction is a 
constant. A perturbation corresponding to such a constant mode generally represents 
a variation of the parameters of the background solution such as A, M, and Q. Thus, 
this mode is relevant to arguments on perturbative uniqueness of a background 
solution. Note that when fc^ = is contained in the full spectrum but does not 
belong to the point spectrum, as in the case = M", it can be ignored without 
loss of generality. 

For modes with k'^ > 0, we can use the vector fields and the symmetric trace-free 
tensor fields defined by 

Si = -^AS, (3-2a) 



1 ^ A „ 1 



k^^ 



DiDjS + -7ijS; S\ = (3-2b) 



n 



to expand vector and symmetric trace-free tensor fields, respectively. Note that Sj 
is also an eigenmode of the operator D ■ D, i.e.. 



while Sij satisfies 



[D-D + k'^-{n- l)K]Si = 0, (3-3) 



(Ai - fc2)§. . = 0, (3-4) 

where A^ is the Lichnerowicz operator defined by 

ALhij := -D ■ Dhij - 2Rikjih^^ + 2(n - l)Khij. (3-5) 

When is a constant curvature space, this operator is related to the Laplace- 
Beltrami operator by 

Kl = -A + 2nK, (3-6) 

In the case of scalar harmonics, the modes with k"^ = uK are exceptional. Given 
our assumption, these modes exist only for = 1. Because is compact and 
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closed in this case, from the identity 

n — 1 — uK 
n k 



Dj&i = — Si , (3-7) 



we have DjS^i = 0. From this, it fohows that J (r'Zy/jE>ij*§>^^ = 0. Hence, E>ij 
vanishes identically. 

We can further show that the second smallest eigenvalue for —A is equal to or 
greater than uK when K > 0. To see this, let us define Qij by 

Qij := LijY = DiDjY - ^lij^Y . 
Then, we have the identity 

QijQ'^ = D\&YDiDjY-YDiAY-RijD^Y)+Y A(A + (n - 1)K) Y-^{AYf . 

For y = S, integrating this identity, we obtain the constraint on the second eigenvalue 

k'^>nK. (3-8) 

For Jf^" = 5"*, the equality holds for the second smallest eigenvalue as we see soon. 

3.1.2. Harmonic functions on S'^ 

Explicit expressions for the harmonic functions on higher-dimensional spaces 
sometimes become necessary to investigate the global structure of perturbations. The 
multiplicity of eigenvalues also has a crucial importance in applying the perturbation 
theory to black hole evaporation. Here, we give such information for harmonic 
functions on 5". 

There are several ways to express higher-dimensional spherical harmonic func- 
tions. For example, we can get expressions in terms of special functions by solving 
the recurrence relation obtained by dimensional reduction. In this paper, we give 
a different approach in which harmonic functions are expressed in terms of the ho- 
mogeneous cartesian coordinates for the Euclidean space containing the unit 
sphere 5". 

Let us denote the homogeneous cartesian coordinates of S"' by il^ = 1, • • • , n+ 
1); ■ il = 1. Then, we can show that satisfies 

DiDjQ"^ = --fijf2^, (3-9a) 

AA^?^ = -Dif2^, (3-9b) 
biQ^D'Q^ = 5^^ - Q'^Q^. (3-9c) 

From these formulae, the following theorem holds. 

Theorem 3.1 (Scalar harmonics on S"). Let us define the function Ya on by 

Ya{n) = aA^-A,n^^ ■■■n^' (3-10) 

in terms of a constant tensor a = (aA^-Ai) (^i, • • • , ^£ = 1, • " " > ^ + !)• Then, the 
following statements hold: 
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1) Ya is a harmonic function on S"" with the eigenvalue 

k"^ = i{£ + n-l), £ = 0, 1, 2, • • • , (3-11) 
if and only if a satisfies the conditions 

aAi-Ai= a^Ai-Ae), (3-12a) 
aA,...A,_,^B = {£>2). (3-12b) 

2) The harmonic functions {Ya} form a complete basis in L^(S'^). Two harmonic 
functions with different values of i are orthogonal and those with the same i 
have the inner product 

(y„,y„0 := J d'^QYaYa' = C{n,lW-''aj,...j,. (3-13a) 

n+l 

3) The multiplicity of the ^-eigenvalue is given by 

^,^^,.,^ (n + 2.-l)(„ + .-2)l ^ (3.14) 

Proof. The first statement follows from 

A„y„ = -InYa + ^ aAi-A.^^"^^ • • • Dk^'^'' ■ ■ ■ D''f2^'' ■ ■ ■ Q^' 



= -m + n- l)y„ + ^Jl^aA,...Ae.,''Bf2^' ■ ■ ■ f^^'-' 



(3-15) 



which can be easily verified with the helps of the identities, Eqs. (3-9). 

Next, all polynomials in the cartesian coordinates for are dense in the 

space of continuous functions in the unit cube in E"^~^^, hence its restriction on S"" 
is also dense in the space of continuous functions on S*". This implies that all the 
harmonic functions of the type Ya are dense in the function space L^(S'"). This 
proves the completeness. The inner product of these harmonic functions can be 
calculated as follows. First, by differentiating the function 

nr') := J d-Qe-^ = Q^., j\e sW^~' Oe^-^^ = ^nj^j0=^J^{r) 

(3-16) 

repeatedly with respect to x"^ (r = {x^)), we obtain 



/ 



(Pf2Yaif2)e''-^ = (-l)^2^x^i • • • x^'aA,...AeF^^\r^)- (3-17) 



Differentiating this £-th times with respect to x^ and putting a;"^ = yield Eq. (3-13). 
The multiplicity formula of the eigenvalue can be easily obtained by just counting 
the linearly independent solutions to Eq. (3-12). □ 



'Perturbations and Stability of Static Black Holes in Higher Dimensions' 15 



3.2. Vector harmonics 

3.2.1. Definition and properties 

A harmonic vector is defined as a vector field on 5" satisfying 

{D ■ i) + kl)Yi = 0; &Yi = 0. (3-18) 

Prom this we can define a symmetric trace-free tensor of rank 2 by 

V^. = -M.V,-), (3-19) 

where the factor 1/ky is just a convention (see below for the case ky = 0). This 
tensor is an eigentensor of the Lichnerowicz operator, 

AlV.j = + {n- 1)K] , (3-20) 

but is not an eigentensor of the Laplacian in general when is not a constant 
curvature space. 

In this paper, we assume that the Laplacian —D D is extended to a non- negative 
self- adjoint operator in the L^-space of divergence- free vector fields on c^", in order 
to guarantee the completeness of the vector harmonics. Because —D-D is symmetric 
and non-negative in the space consisting of smooth divergence-free vector fields with 
compact support, it always possesses a Priedrichs extension that has the desired 
property. With this assumption, k"^ is non-negative. 

One subtlety that arises in this harmonic expansion concerns the zero modes of 
the Laplacian. If =J^" is closed, from the integration of the identity D^{V^ DiVj) — 
D'VWiVj = VW ■ DVj, it follows that DiYj = for fc^ = 0. Hence, we cannot 
construct a harmonic tensor from such a vector harmonic. We obtain the same result 
even in the case in which J^T"^ is open if we require that Y^D^Yj fall off sufficiently 
rapidly at infinity. In the present paper, we assume that this fall-off condition is 
satisfied. From the identity D^DiVj = DiD^Vj -|- (n — l)KVi, such a zero mode 
exists only in the case K = Q. Wc can further show that vector fields satisfying 
DiVj = exist if and only if is a product of a locally fiat space and an Einstein 
manifold with vanishing Ricci tensor. 

More generally, Yij vanishes if Vj is a Killing vector. In this case, from the 
relation 

2k,bjYh = [kl -in- 1)K] Yi , (3-21) 

k^ takes the special value k^ = {n — l)K. Because fc^ > 0, this occurs only for 
K = OT K = 1. In the case K = 0, this mode corresponds to the zero mode 
discussed above. In the case K = 1, since we are assuming that J^" is complete, 
is compact and closed, as known from Myers' theorem, and we can show 
the converse, i.e. that if ky = {n — l)K, then Yij vanishes, by integrating the 
identity = YW^Yij = D^{Y'Yij) + k^Y'^Yij over Furthermore, using the 

same identities 



2D^Vj]D^'V^^ = 2Di{VjD^'V^^) + Vj [-A + (n - 1)K] , (3-22a) 
2L>(iT/j)L>(V^) = 2Di{VjD^W^^) + Vj [-A - (n - 1)K] . (3-22b) 



16 



A. Ishibashi and H. Kodama 



we can show that there is no eigenvalue in the range 0<k^<{n- 1)\K\ where the 
second equahty holds only for if = —1. 

3.2.2. Harmonic vectors on 

We can give explicit expressions for vector harmonics on 5" in terms of the 
homogeneous coordinates using Theorem 3.1. 

Theorem 3.2 (Harmonic vectors on S'"). Let us define the vector field VJJ by 

= bA,...A,;BO^' ■ ■ ■ n^'b'Q'' . (3-23) 

in terms of a constant tensor h = {aAi-At,B){'^ii ' -B = 1, • • • , ra + 1). Then, 
the following statements hold: 

1) is a divergence-free harmonic vector on S'^ with eigenvalue 

kl = ^{^ + n-l)-l, £ = l,2,---, (3-24) 

if and only if the constant tensor b satisfies the conditions 

bAi-Ae;B = , (3-25a) 

bA^-A,_J;B = 0, (3-25b) 
b{A,-Ae;Ae+,) = . (3-25c) 

2) All harmonic vectors { V^^ } form a complete basis in the space of divergence- 
free vector fields on 5". Two harmonic vectors with different values of i are 
orthogonal and those with the same i have the inner product 

{Vb, HO := J dI'f2{Vb)iVi, = C{n,e)b ■ b' , (3-26) 

where C{n,i) is the number given in Eq. (3-13)D 

3) The multiplicity of the ^-th eigenvalue is 

^^^^ > - (^ + l)(^-l)!(n-2)! • 

Proof. For a harmonic vector with the eigenvalue , let us define a set of functions 
on 5*" by = V^DiQ^. Then, from Eq. (3-9) wc obtain Kv^ = -^(kl + l)V^. This 
implies that is a harmonic function with the eigenvalue k^ = ky + l = £{£+n — l). 
Therefore, Vi = VsDiQ^ can be written in terms of a constant tensor 
satisfying the conditions, Eqs. (3-25a) and (3-25b) as in Eq. (3-23). The divergence 
of this expression can be written 

biV' = [-(n + e)bA,...A,;A,^, - eaA,...A,_MAe+,] , (3-28) 

where 

OAi-A^.i := bAi-Ae-ii-^ ■ (3-29) 

Prom this, it follows that the divergence free condition for can be expressed as 

{n + £)b(^Ai-Ar,Ae+i) = ^(^{Ai-Ae-i^AeAe+i) ■ 
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After a short caculation, we find that this condition is equivalent to Eq. (3-25c). 
This proves the first statement. 

The completeness in the second statement immediately follows from the com- 
pleteness of the harmonic vectors. The formula for the inner product can be derived 
by the same method as used for the harmonic scalars. 

Finally, we can show that the conditions, Eqs. (3-25), on b are reduced to the 
conditions on 6^'^...^/.^/ and 6n+iA;-4_i;B' with A[, ■■■ ,A'^,B' = 1,--- ,n, 

b{A[-A',:B') = , (3-30a) 

{£ - 2)bn+l^ B'{A\-A'^_^;A'^_^) + '^K+IA\-A[_.j^^ ;B' , (3-30b) 

and all the other components are uniquely determined from these components. By 
counting the number of linealy independent solutions to these conditions, we obtain 
the multiplicity in the theorem. □ 

3.3. Tensor harmonics 

3.3.1. Definition and properties 

In the Einstein space no condition is imposed directly on the Riemann 

tensor BJ'jki itself, though the Ricci tensor is assumed to be proportional to the 
metric. Hence, in general the Lichnerowicz operator appears in the tensor-type 
perturbation equation instead of a simple sum of the Laplace-Beltrami operator 
and scalars as ^"-dependent part. Thus, we have to use the eigentensors to the 
Lichnerowicz operators to expand tensor-type perturbations on a generic Einstein 
space 

AiTij = \LTij , (3-31) 
Ti = 0, D^ij = . (3-32) 

Note that the Lichnerowicz operator preserves the trace-free and transverse condi- 
tions. When J^f"^ is a constant curvature space, Tjj becomes a harmonic tensor on 

(A + k\)Ti^ = , = Xl- 2nK . (3-33) 

Note also that an Einstein space with dimension equal to or smaller than 3 is always 
a constant curvature space. Further, there exist no symmetric harmonic tensor with 
rank 2 on S'^ and very special ones on and /F. To be precise, the following 
theorem holds: 

Theorem 3.3 (Harmonic tensors on 2-dimensional constant curvature space). Let 
be a two-dimensional closed surface with a constant curvature K. Then, a sym- 
metric harmonic tensor Tjj with rank 2 represents a moduli deformation of and 
exists only when < 0. For {K = 0), Tij is a constant trace-free tensor with 
k"^ = in the chart in which ds'^ = dx^ + dy^, while k"^ = -2K for H^/F {K < 0). 

Very little is known about the spectrum of the Lichnerowicz operator on a generic 
Einstein space. However, we can easily show that 



k^ >n\K 



(3-34) 
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for tensor harmonics on a constant curvature space with the helps of the identities 

2D[iTj]kD^'T^^'' = 2D\TjkD^'T^^^) + Tjk{-^ + nK)T^^ , (3-35a) 

2D^Tj-jkD^'T^^'' = 2D\TjkD^'T^^^) + Tjk{-A - nK)T^^ . (3-35b) 

3.3.2. Harmonic tensors on 5" 

We can give expUcit expressions for the symmetric harmonic tensors of rank 2 
on S*" in terms of homogeneous coordiantcs as in the case of harmonic vectors. 

Theorem 3.4 (Harmonic tensors on S*"). Let us define a tensor of rank 2 by 

Tcij = CA,...At,B,B,^2^' ■ ■ ■ Q^'biQ^^bjQ^^ (3-36) 

where c = {ca^---Ai:B\B2) (^Ij ' ' ' ; ^ii B2 = 1, - ■ ■ , n + 1) is a constant tensor on 
^n+i^ Then, the following statements hold: 

1) Tcij is a harmonic tensor with the eigenvalue 

A;2 = + n - 1) - 2 = 2, • • • ) (3-37) 
if and only if c satisfies the following conditions: 

CAi-At,BiB2 = C(Ai-A(y,BiB2 > (3-38a) 
CAi-Ai.2A^;BiB2=^-, (3-38b) 
C(Ai-Ae;At+i) = C(^Ai-At; A^+i) = . (3-38c) 

2) The set of all harmonic tensors of the form Tdj forms a complete basis for the 

space of trace- free and divergence- free tensors of rank 2 on S"". Two such 
harmonic tensors with different values of i are orthogonal and those with the 
same £ have the innter product 

(Tc, Te/) := J d'^nfcijTc''^ = C{n, i)c ■ c' , (3-39) 

where C{n,i) is the same constant as that in Eq. (3-13)D 

3) For ^ > 2, n > 2, the multiplicity of the £-th eigenvalue for symmetric harmonic 
tensors of rank 2 is 

t c^n^ _ (n + l)(n - 2)(n + e){n + 2£ - l){n + £ - 3)\ 
^^^^ ' - 2^(£-2)!(n-l)! • 

Proof. The first two statements can be easily proved by methods similar to those for 
the harmonic vector. To prove the last statement, lengthy calculations or sophisti- 
cated considerations based on group representation theory are required. See Ref. 23) 
for details. □ 

§4. Tensor-type perturbations 

4.1. Generic background 

We first consider tensor perturbations in an (m + n)-dimensional generic back- 
ground metric, Eq. (2-1), and the energy-momentum tensor, Eq. (2-5). We have 
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already seen in the previous section that hj, ij, 5T^'ij are by themselves gauge- 
invariant. We expand these two in terms of the eigentensors Tjj of the Lichnerowicz 
operator A^, introduced in the previous section as follows: 



.(2). 



h 



(2) 
T ^i 



(2) 

T 2J 



(tt + 2PHt) T 



(4-1) 



The Einstein equations in terms of the gauge-invariant coefficients Ht and tt are 
obtained from Eq. (23) in Ref. 24) with the replacement A;^ — t- Aj, — 2nK, where 
is the eigenvalue of the Lichnerowicz operator. The result is expressed by the single 
equation 

XL-2{n-l)K 



n 



UHt + -Dr ■ DHt 



(4-2) 



We emphasise that this equation holds whenever the background metric is given 
in the form of Eq. (2-1), irrespective to the dimension of ,yV"^, and therefore applies 
to a more general background than that of a static black hole. For example, it has 
been applied to the stability analysis^^^ of Myers-Perry black holes with a single 
rotation, whose metric takes the warped product form of Eq. (2-1), with J^"' being 
the {d — 4)-dimensional unit sphere. 

We also note that Eq. (4-2) is precisely the same form as the equation of motion 
for a massless test scalar field on the same background spacetime if Ht is viewed as 
the scalar field with the angular momentum number being A^ — 2(n — 1)K. 

4.2. Static black hole background 

We turn to the black hole background, Eq. (2-14), for which m = 2. If one 
introduces the master variable # by 



Eq. (4-2) can be put into the canonical form 



(4-3) 



where 



Vt 



AL-2(n-l)i^ + 



nrf n{n — 2)f 



2 4 

In particular, for /(r) given by Eq. (2-15), Vt is expressed as 



(4-4) 



(4-5) 



Vt = 



f 



- lOn + 8 n(n + 2) , , 
Xl H : K — -Xr^ + 



v?M n(3n - 2)Q'^ 



2r 



n— 1 



(4-6) 



Note that since an electromagnetic field is described by a vector field, it does 
not have any tensor-type component. The electromagnetic field enter the equations 
for a tensor perturbation only through their effect on the background geometry. 
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§5. Vector-type perturbations 

5.1. General background case 

Wc have already introduced a basis of vector-type gauge-invariant variables 
{Fj^]\ rj^^^ T-j^ ) in generic background of (m + n)-dimcnsions. We expand these 
gauge invariants in terms of vector harmonics Vj and write the perturbed Einstein 
equations for the expansion coefHcients. 

For generic modes my ■= — {n— 1)K ^ 0, 

F^^^rFaYi, T«=rT„V,, = rVV,,- , (5-1) 

and the Einstein equations reduce to 



1 



r 



^F, = -2kV„ , (5-2) 



^Da{r--'F-) = -K\T. (5-3) 



For the exceptional mode mv = 0, only the following combination 

is gauge-invariant. For this mode we have only a single equation 

_l_^^>(,n+i^(i)^^_2^2^^_ (5-5) 

5.2. Static black hole background case 

In the black hole case with m = 2, using the 2-dimensional Levi-Civita tensor 
€ab, we can rewrite the perturbed Einstein equations as 

Da (r'^+^F^i)) - rnvr^'-'eabF' = -2A^V"+^ea6T^ (5-6a) 

A;^L>a(r"-^F«) = -kV"tt, (5-6b) 

where 

= e'^^Da , (5-7) 

and Fa is defined by the first of Eq. (5-1) with — t- h^^i ■ This should be supple- 
mented by the perturbation of the energy-momentum conservation law 

i^a(r"+V") + ^rVr = 0. (5-8) 

Note that, for my = 0, the perturbation variables hip\ and dT^\ — hence Ht 
and tt — do not exist. The matter variable Ta is still gauge-invariant, but concerning 
the metric variables, only the combination F^^") defined in Eq. (5-7) is gauge invariant. 
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In this case, the Einstein equations are reduced to the single equation (5-6a), and 
the energy-momentum conservation law is given by Eq. (5-8) without the tt term. 

These gauge-invariant perturbation equations can be reduced to a single wave 
equation with a source in the 2-dimcnsional spacetime First, for the generic 

modes my 7^ 0, from Eqs. (5-8) and (5-6b) we find that can be written in terms 
of a variable f2 as 

e^'^Dbf? = r^'-^F" - — r"+V" . (5-9) 

my 

Inserting this expression into Eq. (5-6a), we obtain the master equation 

r^Da {\d<^^\ - = -^r^e-^Dairn) . (5-10) 

For the special modes my = 0, it follows from Eq. (5-8) with tt = that can 
be expressed in terms of a function r(i) as 

r"+Va = e„6L>VW. (5-11) 

Inserting this expression into Eq. (5-6a) with e'^Dc{Fd/r) replaced by F^^^ /r, we 
obtain 

Da{r'''+^F^^y) = -2K^DaT^^\ (5-12) 

Taking into account of the freedom of adding a constant in the definition of t^'^\ we 
have the general solution 

Hence, there exists no dynamical freedom in these special modes. In particular, in 
the source-free case in which r^^^ is a constant and K = 1, this solution corresponds 
to adding a small rotation to the background solution. 

5.3. Static black hole in Einstein- Maxwell system 
5.3.1. Perturbation of electromagnetic fields 

For vector- type perturbations, we note that 5^mn is gauge-invariant as well 
as invariant under a coordinate gauge transformation. As shown in Eq. (2-55), a 
vector perturbation of the Maxwell field can be expressed, in terms of the single 
gauge-invariant variable ^ defined by dss/l^"* = s^Ni on jY"^ ^ as 

b^ah = , 5^ai = Da^Ni , d^ij = ^ (aVj - DjYi^ . (5-14) 
We also expand the current Jj as 

Ji = JYi. (5-15) 

Then we obtain from the second of Eq. (2-52b) the gauge-invariant form for the 
Maxwell equation, 

J—DJr^-'^D^j^) - M±i!iZLl)^^ = - J + rEoPW . (5-16) 
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In order to complete the formulation of the basic perturbation equations, we 
must separate the contribution of the electromagnetic field to the source term in the 
Einstein equation (5-10). The contributions of the electromagnetic field to Ta and 
tt are given in terms of J^/ by 

ri-) = -:^e,,Z)V, 4"'^=0. (5-17) 

Hence, the Einstein equations for the Einstein-Maxwell system can be obtained by 
replacing in Eq. (5-10) by 

r„ = T(^"")+t„, (5-18) 

where the second term represents the contribution from matter other than the elec- 
tromagnetic field. 

5.3.2. Master equations 

Now, generic modes: my 7^ 0, we introduce new master variables by 

.2, 



^± := a±r-"/2 ( [} - + 5±r"/2- V , (5-19) 

V rnv J 



with 



'°--^->-( S[(JVf^i )' '^■^"^^ 

where is a positive constant satisfying 

Z\2 = (n^ - ifM^ + 2n(n - l)mvQ^ ■ (5-21) 

Then, from Eqs. (5-10) and (5-16), we obtain the two decoupled wave equations as 
our master equations: 

Vv+ 

- -Y^± = Sv± , (5-22) 



where 



■^2 ^ (n^ - 2n + 4)K _ njn - 2) ^ n(5n - 2)Q^ ^ /x± 



4 4 4j,2r!.-2 j.n—1 



(5-23) 



and 



fi± = ^M±A, (5-24) 

9 2 n/2 f 

Sv± = -a± ^ e-'Dairn) - b^r^/^'^fJ. (5-25) 

my 

For n = 2, K = 1 and A = 0, the variables #+ and #_ are proportional to the 
variables for the axial modes, z[~^ and ^2"-' given in Ref. 26), and Vy+ and Vy_ 
coincide with the corresponding potentials, V^~^ and V2~'\ respectively. 
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(5-27) 



Here, note that in the hmit Q 0, <P+ becomes proportional to ^ and ^_ to i7. 
Hence, (P+ and represent the electromagnetic mode and the gravitational mode, 
respectively. In particular, in the limit Q — t- 0, the equation for coincides with 
the master equation for a vector perturbation on a neutral black hole background 
derived in Paper I. 

As for the exceptional modes, my = 0, from the definition, Eq. (5-11), of r^^^^ 
and Eq. (5-17), we can express t^^^ as 

r« = -5^ + rW. (5-26) 

Hence, Eq. (5-13) can be rewritten as 

^(1) ^ 2/.^(9^-fW) 

Inserting this into Eq. (5-16), we obtain 

^^,D.(r-^D'^) - J, (2(„ - 1)A- + 5l<=^) ^=-J- « . (5.28) 
Therefore, only the electromagnetic perturbation is dynamical. 

§6. Scalar- type perturbations 

6.1. General background case 

In terms of scalar harmonics S, we expand the scalar- type perturbation variables 

as 

hab = fab^, /ia = -^/aS, hL=2r''HL^, hT = 2r''^8, (6-la) 

STab = TabS , ^Ta = -^ (Pfa + Ta)§ , (6- lb) 

STl = r'^{2HLP + (5P)S , 6Tt = -j^i^HrP + tt)S , (6-lc) 

and Xa = XaS. Here P and all expansion coefficients, e.g. SP, are tensor fields on 
the; -dimensional spacetime 

The basis of the scalar-type gauge-invariant variables introduced in the previous 
section are expanded as: 

= FabS , = 2r^FS , (6-2a) 

= EabS , 17(0),, = rEaSi , (6-2b) 

i:(0) = r^SS , n^o)-. = r^TT^ij . (6-2c) 

The expansion coefficients here Fab, F, ^abj ^ j'T't ^ well as Xa are precisely the 
same as those given in Ref. 24). 

For the exceptional modes with k^ = n ioi K = 1, hr and 6Tt (equivalently 
Ht and tt) are not defined, because a second-rank symmetric tensor cannot be 



24 



A. Ishibashi and H. Kodama 



constructed from S for these modes. In this case, we define F, Fab, ^ab, and El 
by setting X^, = in the above definitions. These quantities defined in this way are, 
however, gauge dependent. These exceptional modes are treated in Appendix C of 
Paper III. 

6.2. Static black hole background 
6.2.1. Maxwell equations 

From now on we consider the static black hole background with m = 2. We 
often use Cab- We expand S''^^^ and S'a'^^ defined in Eqs. (2-57), (2-58) as 

<fS = <f(°), (4S = -^^„(o). (6-3a) 
Then, the Maxwell equations (2-52) are written 

\Da{r^<^) + -^a - ^DaiF^ - 2nF) = eabJ\ (6-4a) 
e«^r»a(r"-Vf,) = -r"-^J, (6-4b) 

with J defined by Ji = rJE>i. Note that from the first of Eqs. (2-51) we have the 
relation 

^ = -^D,{r^'^). (6-5) 

Note also that Eqs. (6-4a) and (6-4b) give the current conservation law 

Dc{r'^J'') = -kr'^-^J. (6-6) 

We find that the gauge invariant variables, Sa and can be expressed in terms 
of the single master variable as 



« = ip-"^ + ' ^"-^ = -^'-^ + - 2nF) , (6-7) 

where Ja has been defined by J" = k'^r~^e"'^ Ji,. By inserting these expressions into 
Eq. (6-5), we obtain the wave equation for 

r-'Da (^) - = -r---Da - ^(F,^ - 2nF) . (6-8) 

The contribution of the electromagnetic field to the perturbation of the energy- 
momentum tensor, Eab, ^a and Sl, are 

2 

Fab , (6-9a) 



yn(em) _ 


^ j.2n ^ 
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nq 
^2n 


F^gab 


v-.(em) _ 


qk 


(Da^ 


' + Ja) 




^(em) ^ 


qk"^ ^ 
j.2n 


2 

nq 

j.2n 


F. 



2j.2n 

(6-9b) 
(6-9c) 
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6.2.2. Master equations 

For generic modes of scalar perturbations, the Einstein equations consist of four 
sets of equations of the forms 

Eab = n'Sab , Ea = K^Ea , El = , Et = k'tt ■ (6-10) 

For the definitions of Eab, Ea-, El and Et, see eqs. (63)-(66) in Ref. 24). Introducing 
the perturbation of the energy-momentum tensor as 

n— 1 2 

Sab = r--'^\Kb - Sa = '-^{^a - 4^^")), Sl = r--'n\EL - E^r^) , 

(6-11) 

and 

St = (6-12) 
the conservation law for the energy-momentum tensor is given by the two equations 

lW)-g.+ <"-^4*;-"^' 5. = 0, (6.13a) 



^A(r-S;) + !Ls, - n^S, = *^^J.. (6.13b) 

Now we introduce X, Y and Z defined by 

X = r'^-'^{Fl -2F), Y = r"-2(F; -2F), Z = r^-^F^, (6-14) 

as in Paper I. After the Fourier transformation with respect to the Killing time coor- 
dinate, t, of the black hole background, Eq. (2-14), the perturbation equations above 
can be reduced to a system consisting of three first order linear differential equations 
for X, Y, Z and a single linear algebraic constraint on them, with inhomogeneous 
terms given by £/, Ja, Sab, Sa, St- As performed in Paper I, this constrained system 
can be further simplified to a single second-order ODE for a scalar field <P given by 



with a source term, where 



, n(n+l)M n^Q^ 
H = '^ + ^:^i (6-16) 

m = k^-nK. (6-17) 

Thus obtained ODE for ^ and the Maxwell equation for s^^ , given by Eq. (6-8), 
form a coupled second-order ODEs with source terms. Then, by introducing new 
master variables, given as a linear combination of ^ and ^ below, we obtain 
the following two decoupled master equations for scalar-type perturbations: 
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where 

#± := a±^ + 6±^/ , (6-19) 

with 

K.M^(^Q. '""^:i^^"'" Q. '"";:itt'"^ ). (6.20a, 

(a_,6_)=((„+l)(A/ + W-^.-^) , (6-20b) 



and where /x is a positive constant satisfying 

f^' = M- + -^^. (6-21) 
If we define the parameter S by 

H = {1 + 2mS)M , (6-22) 
the effective potentials Vs± are expressed as 



with 



^^^^_ n(n + l) ^^^ ^_ =^ + !^(!^_ti)(l + ^<5)^^ (6.24a) 



and 



a;2 



U+ = [-4n^(n + 2)(n + l)^^^^^ - 48n^(n + l)(n - 2)fe 

— 16(n — 2)(n — 4)] y — (3n-2)(n + 1)^(1 + m(5)x^ 
+4(^^n2(n + 1)^ {(n + l)(3n - 2)mS + 4n^ + n - 2} 
+46{n + 1) {(n - 2)(n - 4)(n + l)(m + n^/i:)5 - 7n'' + Tn^ - 14n + 8} 
+ {I6(n + 1) (-4m + 3n2(n - 2)^:) 5 - 16(3n - 2)(n - 2)} x 
+64m + 16n(n + 2)K , (6-25a) 

U- = [-4n^{n + 2)(n + 1)^(1 + mSfx^ + 48n^(n + l)(n - 2)m(l + m(5)x 
-16(n - 2)(n - A)rr?] y - n^{3n -2){n + lfS{l + mSfx"^ 
-iri^in + 1)2(1 + mSf { (n + l)(3n - 2)m5 - n^} x^ 
+4(n + 1)(1 + mS) {m(n - 2)(n - 4)(n + l)(m + n^i^)*^ 



X 



2 



+4n(2n2 - 3n + 4)m + n2(n - 2)(n - 4)(n + l)if } 
-16m {(n + l)m (-4m + 3n'^{n - 2)K) 6 
+3n(n - 4)m + Sn^{n + l)(n - 2)K} x 

+64m^ + 16n(n + 2)711^ K . (6-25b) 
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The source terms Ss± are given by 

Ss± = a±S^ + b±Sj^ , 

where S^ = 5^1^=0 with 
/ 



(6-26) 



S^ 




+ 2nrfJr + 2k^'^ + 2nf'^ 
iuj iu 



_HSr-^^-2nf'-^-2nrfSr 



H ioj 



IUJ 



Ps3rSl 
H iuj 



+ 2r 



ioj 



+ 2r^S: 



(6-27) 



and 

Sa = 



2n\n-lfzp , 2 

2TJ 



Jt 



IUJ 



fJr \ 2(n - 1)^0 



iujH 



finfSt-rSl 



(6-28) 

Here Psi,Ps2 and P53 are polynomials of x,y and z, whose explicit expressions are 

Pgi = [-4n^z + 2n^(n + l)x - An{n - 2)m] y 

+ {2n^(n - l)x + 4n{n - 2)m + 4n^(n - 2)K} z 
-n^iji^ - l)x^ + {-4n(n - 2)m + 2n^(n + l)K] x 
+4m^ + An^mK , 

[Gn'^z — n^(n + l)(n + 2)x + 2n{n — 4)m] y 
-2n'^z'^ + {n^{3n^ - n + 2)x - 4n(n - 2)m - 6n^(n - l)K} z 



(6-29) 



'{n + l).'r2 + {n(3n - 7)m + ra^(n^ - l)K} : 



— 2m — 2n(n — l)m-fC , 
P53 = -2n^{3n - 2)z + n^{n + \)x - 2{n - 2)m . 



(6-30) 
(6-31) 



Note that the master variable coincides with that for the neutral and source-free 
case in Paper I for Q = and St = 0. Also note that the following relations hold: 



q2 ^ (n + lfM^5{l + mS), 
H = H+H^ . 



(6-32) 
(6-33) 



Prom these relations, we find that (5 = corresponds to (5 = 0, and in this limit, 
<P_ coincides with and its equation coincides with the mater equation for the 
master variable # derived in Paper I. Hence, and #-|- represent the gravitational 
mode and the electromagnetic mode, respectively. For n = 2, = 1 and A = 0, these 
variables <P+ and #_ are proportional to the variables for the polar modes, and 
Z2^\ appearing in Ref. 26), and V5+ and Vg- coincide with the corresponding 
potentials, y/"*"^ and Vg^^ respectively. 
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For the exceptional modes, the last of Eq. (6-10) is not obtained from the Ein- 
stein equations. However, this equation with tt = can be imposed as a gauge 
condition, as shown in Paper I. Under this gauge condition, all equations derived in 
this subsection hold without change. However, the variables still contain some resid- 
ual gauge freedom. See Appendix D of Paper HI, for how to eliminate the residual 
gauge freedom to extract physical degrees of freedom. 

§7. Lovelock Black Holes 

When Einstein derived the field equations for gravity, he adopted the three 
requirements as the guiding principle in addition to the principle of general relativity. 
The first is the metric ansatz that gravitational field is completely determined by the 
spacetimc metric. The second is that the energy-momentum tensor sources gravity, 
hence it is balanced by a second-rank symmetric tensor constructed from the metric. 
The last is the requirement that the field equations contain the second derivatives of 
the metric at most and are quasi- linear, i.e., the coefficients of the second derivatives 
contain only the metric and its first derivatives. These requirements determine the 
gravitational field equation uniquely. 

In 4-dimensions, we can obtain a similar result even if we loosen the third con- 
dition and require only that the field equations contain second derivatives of the 
metric at most. To be precise, gravity theories satisfying these weaker requirements 
are the Einstein gravity and the f{R) gravity, the latter of which is mathemati- 
cally equivalent to the Einstein gravity coupled with a scalar field with a non-trivial 
potential. 

When we extend general relativity to higher dimensions, however, the difference 
between the two versions of the third requirement becomes important. In fact, 
we require the stronger version, we obtain the same field equations for gravity in 
higher dimensions. In contrast, we require the weaker version, we obtain a larger 
class of theories that contains the higher-dimensional general relativity as a special 
case. Such extensions to higher dimensions were first studied systematically by 
D. Lovelock in 1971.^7) What he found was that the second-rank gravitational 
tensor Emn balancing the energy-momentum tensor Tmn is a sum of polynomials 
in the curvature tensor and that the polynomial of each degree is unique up to a 
proportionality constant. The physical importance of such an extension was later 
recognised when B. Zwiebach^^^ pointed out that the special quadratic combinations 
of the curvature tensor named the Gauss-Bonnet term naturally arises when we 
add quadratic terms of the Ricci curvature to the quadratic term in the Riemann 
curvature tensor obtained as the a' correction to the field equations in the heterotic 
string theory to obtain a ghost-free theory. In this section, we briefly overview 
the Lovelock theory, its static black hole solution and perturbation theory recently 
developed by Takahashi and Soda.^^^'^^^ 

For notational convenience, we introduce an orthonormal basis on and 
throughout this section, we use upper case latin indices in the range A, B, . . . , J 
to label the 1-form 6^, {A = 0, ■ ■ ■ , D — 1), oi the basis, while we use upper case 
latin indices in the range K,L,M,N,... to denote tensors on ^ as in the rest of 
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this chapter. 

7.1. Lovelock theory 

In 1986, B. Zumino^°) pointed out that the class of theories obtained by Lovelock 
has a natural mathematical meaning. That is, the Lovelock equations for gravity 
can be obtained from the action that is a hnear combination of the terms each of 
which corresponds to the Euler form in some even dimensions: 

. [P-l)/2] 
A;=0 

= h^D , (7-2) 

where is the curvature form with respect to the orthonormal 1-form basis 6^ 

(A, S = 0, • • • , D - 1), e"^-^ = 6"^ A--- A9^, f2D is the volume form, and 

4 = ^6^i::illRA,A,^'''' ■ ■ ■ Ra,,_,a,,^''-'^'' . (7-3) 
The explicit expressions for small values of k are 

/o = 1 , (7-4a) 
h = R. (7-4b) 
l2 = R^- 4RiR^ + Rab^'^Rcd^'' , (7-4c) 
h = R^- 3R{-4RiRi + Rab^^Rcd^^) + 24i?^i?gi?5 + 3Rab^^ Rcd^'' Ref^^ ■ 

(7-4d) 

Hence, if we require that the theory has the Einstein theory in the low energy limit, 
ao and a\ are related to the cosmological constant A and the Newton constant 
= 8ttG as 

"0 = ai = (7-5) 

Further, if the Gauss- Bonnet term ^2 comes from the 0{a') correction in the het- 
erotic string theory, 02 > 0. 

From the Bianchi identity &Mab = 0, the variation of the Lagrangian density 
can be written 

= '^^*) + P-2^fc-l)! ^^^--^^^^^-^"-^^^^^'^' A • • • A e^--^n-., 

= d{.) + [kTC^^sAS^^'^'nA + {-lf-'E(%deZe^} Qu , (7- 

where oj^b is the connection form with respect to 9^, ca is the vector basis dual to 
6^, S) \s the corresponding covariant exterior derivative, = is the torsion 
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2-form, and 

E^'^B = -YAt-tl^MM^''"' ■ ■ ■ . (7-7b) 

Hence, if we treat 0^ and lo^b as independent dynamical variables, the field equa- 
tions are given by 

^«fe£;('=):^ = 0, (7-8a) 

k 

J2kakT(>'^ic = ^- (7-8b) 
k 

If we require the connection to be Riemannian, the second equation becomes trivial 
due to the torsion free condition T^q = 0. Examples of the explicit expressions for 

E^''^^ and T^'^^ic are 

E('^ = -S^ , (7.9a) 
E^^^ = 2Ri - R5^ , (7-9b) 

+SR^R% + 4i?^^i^2^3^BCiC2C3 , (7-9c) 
T^Hc = T^C + ^SfBT8]D , (7-9d) 
T'^'^ic = 86f^{-2RE]TD + R''''''c]dA. + ^^c] " 2i?gr^^^,) 

+12(i?^^BcTD - 2Rfj,Tc^ - 2R^''^d,[bT^{d, ' ^'^RdTEc) ■ (7-9e) 

7.2. Static black hole solution 

7.2.1. Constant curvature spacetimes 

In general relativity, a constant curvature spacetime is always a vacuum solution 
and the curvature is uniquely determined by the value of the cosmological constant. 
This feature is not shared by the Lovelock theory. In fact, the Lovelock theory does 
not allow a vacuum constant curvature solution for some range of the coupling con- 
stants {ofc}, and have multiple constant curvature solutions with different curvatures 
for other range of the coupling constants. 

To see this, let us insert the constant Riemann curvature 

Rmnlk = ^{gMLQNK — guKgNL) (7- 10) 
into the field equation (7- 7b). Then, we obtain 

P(A) = 0, (7-11) 

where 

nX):= E (7-12) 



..0 iD-2k-l)\ 
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For general relativity for which = for A: > 2, this equation has a unique solution. 
In contrast, when D > 4 and 7^ (A; > 2), the equation can have no solution or 
multiple solutions depending on the functional shape of P{X). 

7.2.2. Black hole solution 

Now, let us look for spherically symmetric black hole solutions. Because the 
Birkhoff-type theorem holds for the Lovelock theory except for the case in which 
P(A) = has a root with multiplicity higher than one,^^)'^^) we only consider static 
spacetimes whose metric can be put into the form 

ds^ = -f{r)dt' + r'^dul , (7-13) 

where da^ represents the metric of a constant curvature space with sectional cur- 
vature K. For a spherically symmetric solution, K = \. However, because the 
argument in this section holds for any value of if, we consider this slightly general 
spacetime. 

Then, the non-vanishing components of the curvature tensor are given up to 
symmetry by 



2 V / 2/2 ; 4/ ' 

-^OiOi = -^9i3 ' ^liii = ~2^9ij , (7- 14b) 

Rijki = X{gikgji - gugjk) , (7- 14c) 

where X{r) := {K — h{r))/r'^. Inserting these into Eq. (7- 7b), we find that the field 
equations reduce to 

(r"+ip(X(r)))' = 0, p(i)(X(r))(/(r)//i(r))' = 0. (7-15) 

The first of these is integrated to yield 

^W0) = ;^, (7-16) 

where C is an integration constant. This determines the function h{r) implicitly. In 
particular, for C = 0, we have h{r) = K — \r^ for each solution to P{\) = 0. If 
P(^)(A) 7^ 0, after an appropriate scaling of t, /(r) = h{r) follows from the second 
of the above field equations: 

^ _ ^^2)^^2 ^ + r^dal . (7-17) 

This represents a constant curvature spacetime with sectional curvature A irrespec- 
tive of the value of K, as is well known. 

This implies that for C 7^ 0, we have in general multiple solutions corresponding 
to multiple solutions to P(A) = 0. Each solution approaches a constant curvature 
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spacetime with sectional curvature A at large r asymptotically. For these solutions, 
we can always put /(r) = h{r) and the metric can be written^^^ 



(7-18) 



In general, the constant C is proportional to the total mass M of the system and 
positive if M > 0. We can easily show that X{r) changes monotonically with r from 
infinity to some value of r where the metric becomes singular. This singularity may 
or may not be hidden by a horizon depending on the functional shape of P{X). In 
the former case, we obtain a regular black hole solution. 

7.3. Perturbation equations for the static solution 

The linear perturbation of the Lovelock tensor, Eq. (7-7b), in general reads 

[(D-l)/2] 



-6 



k=l 



2k "NNi-N2k 



Inserting Eq. (7-14) with /(r) = h{r) = K — X{r)r'^ into this yields^^^ 

vT' 



-SRij'^ - 2T5R, 



ir 
ir ) 



„n-l 



„n-l 



„n-l 



„n-l 



n — 1 

SE: = -2T6Rii', 
rT' 
n - 1 

-5Raj'^ +2T5Ra^\ 

n — 1 
2rT' 



6EI 



2r^T" 



6E' 



n 



SRaj""' + 



(n-l)(n-2) 



6Rr'^ 



2rT' 

2TSRtr'' + -5Rak'''' + 

n — 1 



(n-l)(n-2) 



SRkl 



kl 



where 



(7-19) 

(7-20a) 
(7-20b) 
(7-20c) 

(7-20d) 

(7-20e) 

(7-21) 



T{r) :=r"-^pW(X(r)). 

Thus, we can obtain perturbation equations for the metric in the Lovelock theory 
simply by calculating the perturbation of the curvature tensor. 

7.3.1. Tensor perturbations 

For tensor perturbations, the metric perturbation can be expanded in terms of 
the harmonic tensor T.^ as Eq. (4-1). The non- vanishing components of the curvature 
tensor for this type of perturbations read 



SRai^' = - UHt + -Dr ■ DHt 



SRik 



jk 



f 

-{n-2y—H'r + 
r 



2_K+^ ■ 

o 



(7-22a) 
(7-22b) 
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Hence, from the above expression for dEj , we obtain the fohowing wave equation for 
If we introduce the mode function !P'(r) by 



(7-23) 



HT{t,r) 



-iu}t 



r^/T{r) 

this wave equation can be put into the standard form 



dr'i 



with the effective potential 



(7-24) 



(7-25) 



(7-26) 



(n-2)r T' rs/T' drl ' 

where dr* = dr/f{r) as in the Einstein black hole case. 
7.3.2. Vector perturbations 

For vector perturbations, the perturbation of components of the curvature tensor 
that are relevant to the field equations can be expressed in terms of the basic gauge- 
invariant quantities as 



5R 



-aj 



_ (n-2)fc ZJV 



n — 1 



J' 

P — (n 



SRab 



ib 



(n — l)r 



l)K 2{K-f) + rf ■ 

-fa H Ja 



(7-27a) 
(7-27b) 
V\(7-27c) 

(7-27d) 



Inserting these into Eq. (7-20d), we obtain the following equations for the gauge- 
invariant variables: 

lz.'(.=rF«;))+r'i!^^i^F. = 0. (7.28a) 



-L>«(rT'F„) = 0. 



(7-28b) 



The gauge-dependent residuals in the expressions for the curvature tensor cancel 
exactly owing to the identity 



(7-29) 
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obtained from the background equation (r'^'^^P^X))' = 0. We can easily confirm 
that for general relativity for which T = r"~^/(2K^), these reduce to Eqs. (5-2) and 
(5-3) with no source terms. 

A master equation for vector perturbations in the Lovelock theory can be derived 
in the same way as that in general relativity. First, the second perturbation equation 
implies the existence of a potential f2 in which Fa can be expressed as 



(7-30) 



Inserting this into the first perturbation equation, we easily find that it is equivalent 

to 



rTDa 



1 



- {n-l)K 



Q = Q. 



7.3.3. Scalar perturbations 

For scalar perturbations, we have 



(7-31) 
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r 2 



(7-32a) 
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^D^F"^ - —Dr ■ DF" 
r 2r " 



+ + 2^^) ^» - + ^'^^^ (7) } ^)^(^-32b) 

2(^2 _ nK) ^ 



D^rD^r 



Fab--Dr''DaF + 



K-f 



S}S-{n-2)^FS), 



ab 



0.3 



{n-l)k 



F^ + ^DaF- 



Dbv 



2{K^f) + rf 
2r2fc2 



DaHT 



(7-32c) 
(7-32d) 
§\ (7-32e) 
(7-32f) 



|_nF« + D'^D^'Fab + ^F^ + S . 

Inserting these into Eq. (7-20), we obtain 

-^(f;)' - f + — / + ^t) f; + 2nfTF" 

r \ r r r'^ J 

S = , (7-33a) 



+ ( Mr + 2nfT' + nfr] F' - 2T'^^-^^F 
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r^'-^SE: 



Tnf 
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2nfT., nfT ff , T' 
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f 



F + nfT{'j + 2-]F' - 2T 
1 



^,k'^-nK 



Dh -F^ ]+TDa( -Fl; + 2T'DaF 



S = 0,(7-33c) 
= , (7-33d) 



P /T' 
--—i—F^ + 2T"F]8) 
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( f'TY k'^T'\ 

+ F» + D'^D\TFab) 
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(7-33e) 



As was first shown by Takahashi and Soda,-*^^) we can reduce these equations to 
a single master equation in terms of the master variable <P defined by 



= r(# + 2F), 



as 



where 



(7-34) 
(7-35) 



A = 2k^ - 2nf + nrf, 
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f 



nr 
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rik^T + nrfT') (2 



{AT)' T" 



AT 



(7-36a) 

^)-n(rVr')' . (7-36b) 



The other gauge-invariant variables are expressed in terms of <P as 



V.2 A 

F; = -—^^ + 2rF'-—F 
2rT" 

F! = -FI -F. 



(7-37a) 
(7-37b) 
(7-37c) 



§8. Stability Analysis 



8.1. Stability criterion and S-deformation 

With the decoupled master equations in hand, we are ready to study the stability 
of generalized static black holes with charge and cosmological constant in Einstein- 
Maxwell theory, whose metric is given by Eqs. (2-15) and (2-14). We consider only 
the stability in the static region outside the black horizon. This region is represented 
as r > for A < and r// < r < rg for A > 0. Such a region exists only for restricted 
ranges of the parameters M, Q and A. [See Appendix A of Paper III for details.] 
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We have seen before that for any perturbation type, the master equations for 
perturbation in the static region are reduced to an eigenvalue problem of the type 

uj'^^ = A^, (8-1) 

where A is the derivative operator 

dr 

^ = -^ + nr); dr. = -, (8-2) 

with V{r) being equal to Vrir), Vv±{r) or Vs±{r). The operator A is self-adjoint (by 
imposing suitable boundary conditions if necessary) in the standard square integrable 
function space L?'{r^,dr^). Therefore, if its spectrum is non- negative, there is no 
exponentially growing mode among physically acceptable (i.e., normalizable) modes, 
implying the stability of the black hole. 

When A > 0, the static region where A is defined is globally hyperbolic and the 
range of is complete. In this case the operator A is essentially self-adjoint with 
the domain of smooth functions of compact support, denoted hereafter by C^(r*), 
and has the unique self-adjoint extension called the Priedrichs extension Ap, which 
is given by taking the closure of (A, C^(r*)) and is known to have the same lower 
bound of the spectrum of A with domain CQ°(r*). 

However, when A < 0, r* has an upper bound, and whether A becomes essentially 
self-adjoint depends upon the asymptotic behavior of the potential V{r) in A, and 
thus upon the type of perturbations as well as the spacetime dimension. When A 
is not essentially self-adjoint, there are infinitely many different choices of boundary 
conditions that make A self-adjoint, and the spectrum of a self-adjoint extension 
depends upon the associated boundary conditions at the upper-bound of r* (r — >■ oo). 
In particular, even if A with C^(r*) is positive-definite, its extension can admit a 
negative spectrum, depending on the choice of boundary conditions. This could be 
the case for the vector-type of electromagnetic and gravitational perturbations in 
n = 2 and for the scalar-type of electromagnetic and gravitational perturbations 
in n = 2, 3, 4, as shown in a simple, massless case by inspecting the asymptotic 
behavior of A}^^ For these cases, we must specify a boundary condition for <P at 
r = oo, and in the following we simply adopt the setting zero Dirichlet condition, 
— 7- as r — )• OO, which corresponds to the Friedrichs extension, Ap. 

Now suppose (p{r) is a smooth function of compact support contained in r > ru 
(or r// < r < Tc for A > 0). Then, we can rewrite the expectation value of A, 
(#, A<P), as 

{^,A0) = Jdr, ' 



dr* 



+ V\<P\A. (8-3) 



Note that no boundary terms appear as <P G CQ°(r*). When the static region is 
globally hyperbolic, we have the Friedrichs extension Ap- Then, in order to show 
the stability, it is sufficient to show the positivity of the right-hand side of Eq. (8-3) for 
e Co°(r*),'^) since Ap has the same lower bound of (A, C^(r*)). In particular, this 
condition obviously is satisfied if V is non-negative. For 4-dimensional Schwarzschild 
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black hole, this is indeed the case. However, in higher dimensions, V is in general 
not positive definite, and it is far from obvious to know whether A is positive. One 
powerful method to show the positivity of A beyond such a simple situation is the 
procedure called the S- deformation of V in Papers II and III, in which we deform 
the right-hand side of Eq. (8-3) by partial integration in terms of a function S as 

= J dr, (\D0f + V\^^\^ , (8-4) 

where 

D = A. + S, V = V + f^-S\ (8-5) 

and ^ G CQ°(r*). Our task is now to find a suitable function S which makes the 
effective potential V positive. In the following we quote the results of the stabihty 
analysis of Papers II and III. 

One might worry that our boundary conditions ^ — t- in a neighbourhood of 
the horizon r — > r^j would be too strong. However, if the horizon is non-degenerate 
and accordingly admits a bifurcate surface and once the stability is shown under 
our boundary conditions, then we can conclude by applying the theorem of Kay 
and Wald^^-* that the black hole spacctimc is stable under perturbations that are 
non- vanishing at the bifurcate surface. Below we shall examine the stability for each 
type of perturbations. 

8.2. Tensor perturbation 

The stability analysis of higher dimensional static vacuum black holes under 
tensor type perturbations has first been examined by Gibbons and Hartnoll^^^ with 
special interests in the case where J^"" is a generic Einstein manifold. In the following 
we review our analysis in Paper III, which generalize their results. 

Consider the potential, Eq. (4-5). As discussed in Paper II, with the choice 

S = -^, (8-6) 



the iS-deformation yields 



Vt = ^ [\l - 2(n - 1)K] , (8-7) 



irrespective of the r-dependence of /(r). Therefore Vr becomes positive if 

XL>2{n-l)K. (8-8) 

In particular, this immediately guarantees the stability of max;imally symmetric black 
holes for K = 1 and K = 0, since is related to the eigenvalue of the positive 
operator —D^Di as \l = k"^ + 2nK when J(^"' is maximally symmetric. As for 
the case K = —1, Vt could be negative even in the maximally symmetric case, if 
< k^ < 2. This is however not possible, due to the bound, Eq. (3-34). Therefore 
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we conclude that the maximany symmetric black holes considered here are stable 
under tensor type perturbations. 

Note that the condition, Eq. (8-8), is just a sufficient condition for stability, 
and it is not a necessary condition in general. In fact, for tensor type perturbation, 
we can obtain stronger stability conditions directly from the positivity of Vr if we 
restrict the range of parameters. For example, for if = 1 and A = 0, it is easy to see 
that Vt is positive if 

n(n - 1)VM2 -Q2 

Al + 2 - 2n + ^ > (8-9) 

M + - Q2 - ^ ^ 

for M2 > Q2 > 8„(^„ _ i)M2/(3n - 2)2, and 

- lOn + 8 ^ 
Al + ^ >0 (8-10) 

for < 8n(n — l)M2/(3n — 2)2. Thus, if we do not restrict the range of Q^, we 
obtain the same sufficient condition for stability as Eq. (8-8), but for the restricted 
range < 8n{n — l)M'^/{Sn — 2)2, we obtain the stronger sufficient condition, 
Eq. (8- 10), which coincides with the condition obtained in Paper II for the case 
Q = 0. 

Similarly, for K = —1 and A < 0, if we restrict the range of A to 

(n + l)M\^/ (n2-l)M2 



we obtain from Vr > a sufficient condition for stability stronger than Eq. (8-8), 

XL + Sn-2 = k^ + n-2>0. (8-12) 

This condition is sufficient to guarantee the stability of a maximally symmetric black 
hole with K = —1 for n > 2. However, if we extend the range of A to the whole 
allowed range, then Eq. (8-8) becomes the strongest condition that can be obtained 
only from Vr > 0. 



8.3. Vector perturbation 
For the choice of 



S = ^, (8-13) 



V = Vv± in Eq. (5-23) are deformed to 



mv + 



(n2 - 1)M±A 



j.n—1 



(8-14) 



mv = kl-{n-l)K. (8-15) 

It follows from my > that VV+ is positive definite, and therefore static charged 
black holes are stable under the electromagnetic mode of vector type perturbation. 
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n= 




. n= 




n= 


1 // — ■"^ 


2 


^ r/rH 



K=i,'k<0, Q^/M^=0.1, 1=2 



Fig. 1. Examples of Vs+ for K = 1 and A < 0. 



For the gravitational mode, we have 



Vv- 



r2 (n2 - 1)M + A 



, 9 , 2n(n- 1)Q2 



•16) 
•17) 



Since /i is a monotonicahy increasing function of r, Vy- is positive if and only 
if h{rH) > 0. Therefore Vy- may become negative. In the case of A > 0, the 
background spacetime contains a regular black hole only for K = 1, and the static 
region outside the black hole is given hy rfj < r < rc (< +oo). In this region, it 
turns out that h > and therefore the black hole is stable. In the case of A < 0, 
under the condition that the spacetime contains a regular black hole, we have the 
relation 



h > v^(n2 - 1)2M2 + 2n(n - l)mvQ^ - y^{n^ - IfM^ - AKn{n - l)2g2 



•18) 



It follows that for K = 0,1, h > 0. As for K 
inequality could become negative if A;^ < n — 1. 



: —1, the right-hand side of this 
This is however not possible as 



shown just below Eq. (3^22b). Therefore, h > also for K = —1. We conclude that 
the black holes considered here are stable under vector type perturbations. 

8.4. Scalar perturbation 

By applying the S'-deformation to Vs+ with 



S 



dr 



n/2-ljj 



we obtain 



2r'^H. 



[{n-2){n + l)6x + 2] 



•19) 



•20) 
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Fig. 2. Examples of Vs-- 



Since this is positive definite, the electromagnetic mode ^4. is always stable for any 
values of M, Q and A, provided that the spacetime contains a regular black hole, 
although Vs+ has a negative region near the horizon when A < and Q'^/M'^ is small 
(see Fig. 1). 

Using a similar transformation, we can also prove the stability with respect to 
the gravitational mode ^_ for some special cases. For example, the S-deformation 
of Vs- with 

5 = _L^, h^ = r^/^-^H. (8-21) 
n,_ dr 

leads to 

Vs- = iTTTT- [2"^ - (n + l)(n - 2)(1 + m5)x] . (8-22) 
2r^/i_ 

For n = 2, this is positive definite for m > 0. When = 1, A > and n = 3 or 
when A > 0, Q = and the horizon is S"^, we can show that V^^ > 0. Hence, in 
these special cases, the black hole is stable with respect to any type of perturbation. 

However, for the other cases, Vs- is not positive definite for generic values of 
the parameters. The 5-deformation used to prove the stability of neutral black holes 
in Paper H is not effective either. This is because Vs- is negative in the immediate 
vicinity of the horizon for the extremal and near extremal cases, as shown in Fig. 2, 
and the S'-deformation cannot remove this negative region if 5 is a regular function 
at the horizon. Hence, the stability problem for these generic cases with n > 3 is 
left open. 

Our results are summarized in Table I. As shown there, maximally symmetric 
black holes are stable with respect to tensor and vector perturbations over the entire 
parameter range. In contrast, for scalar type perturbation, we were not able to prove 
even the stability of asymptotically flat black holes with charge in generic dimensions, 
due to the existence of a negative region in the effective potential around the horizon 
in the extremal and near extremal cases. 

In this regard, it should be noted that by extensive numerical studies, Konoplya 



'Perturbations and Stability of Static Black Holes in Higher Dimensions ' 41 

Table I. Stabilities of generalized static black holes. In this table, "D" represents the spacetime 
dimension, n + 2. The results for tensor perturbations apply only for maximally symmetric 
black holes, while those for vector and scalar perturbations arc valid for black holes with generic 
Einstein horizons, except in the case with K = 1,Q — 0, X > and D = 6. 





Tensor 


Vector 


Scalar 


Q = 




Q = 




(5 = 


Qt^O 


K = 1 


A = 


OK 


OK 


OK 


OK 


OK 


£» = 4,5 OK 
D>67 


A > 


OK 


OK 


OK 


OK 


D < 6 OK 
D > 7 ? 


D = 4, 5 OK 
D>6? 


A < 


OK 


OK 


OK 


OK 


D = 4 0K 

D>57 


D = 40K 

L> > 5 ? 


K = 


A < 


OK 


OK 


OK 


OK 


D = 4 OK 
D>5? 


D = 40K 
£) > 5 ? 


K = -1 


A < 


OK 


OK 


OK 


OK 


D^40K 
D > 5 ? 


D = 4 OK 
D > 5 ? 



and Zhidenko^^) have shown that Schwarzschild-de Sitter black holes are stable in 
D := 2 + n = 5, ...,11. They also found that charged Schwarzschild-de Sitter black 

holes can be unstable if the electric charge and cosmological constant are large enough 
in D > 7.'^"^^ As for charged asymptotically anti-dc Sitter (AdS) black holes, they 
have found no evidence for instability in D = 5, 11.^^^ 

8.5. Lovelock black holes 

As we have seen in section 7, there are known exact solutions of static black 
holes in Lovelock theory, and the master equations for all types of perturbations 
of static vacuum black holes in general Lovelock theory have recently been de- 
rived by Takahashi and Soda."*^^-* Using the master equations, they have found that 
an asymptotically flat, static Lovelock black hole with small mass is unstable in 
arbitrary higher dimensions; it is unstable with respect to tensor type perturba- 
tions in even-dimensions^^^'"^^) and with respect to scalar type perturbations in odd- 
dimensions.^''^ The stability under vector type perturbations in all dimensions has 
also been shown*^^) by applying the S'-deformation technique. 

In fact, such an instability against tensor and scalar type perturbations, as well 
as the stability under vector type perturbations, have already been indicated by 
earlier work*^-*^)"^^) performed within the framework of second-order Lovelock the- 
ory, often called the Einstein-Gauss-Bonnet theory. For such a restricted class of 
Lovelock theory-though most generic in D = 5,6, the master equations for metric 
perturbations have previously been derived by Dotti and Gleiser.^^^'^^) A numerical 
analysis of the (in)stability of static black holes in Einstein-Gauss-Bonnet theory in 
dimensions D = 5, 11 has been performed in Ref. 45). 

It is interesting to note that the instability found in small Lovelock black holes 
is typically stronger in short distance scales rather than long distance scale/low 
multipoles as one may expect. For example, for tensor type perturbations, there 
appears the eigenvalue of tensor harmonics on the horizon manifold as an overall 
factor in the effective potential term of Eq. (8-4), and it is therefore always possible 
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to make the potential term dominant by taking a sufficiently large eigenvalue of 
tensor harmonics. This implies that if the effective potential term can be negative, 
the right-hand side of Eq. (8-4), as a whole, can also be negative. This is shown to 
be the case when the mass is sufficiently small. A similar argument also applies to 
the case of scalar type perturbations at higher multipole moments. "^^^ 

§9. Summary and Discussions 

We have reviewed a gauge invariant formalism for gravitational and electromag- 
netic perturbations of static charged black holes with cosmological constant in higher 
dimensions and, as an application, the stability analysis using the master equations 
derived in the developed formalism. In section 2, we have started with considering 
a fairly generic class of background spacetimes defined by a warped product of an 
m-dimensional spacetime and an n-dimensional internal space, where the latter cor- 
responds to the horizon cross-section manifold. We have explained how to decompose 
tensor fields in the background spacetime from the viewpoint of the internal space, 
and have seen that second-rank symmetric tensor fields or metric perturbations are 
decomposed into tensor, vector and scalar-types. We have then constructed mani- 
festly gauge invariant variables for each type of perturbations. After that in section 3, 
we have introduced harmonic tensors on the internal space to expand the gauge in- 
variant variables in terms of them. We have presented several theorems concerning 
basic properties of harmonic tensors, and also given explicit expressions of the har- 
monic tensors in terms of homogeneous coordinates. In subsequent sections 4-6, we 
have briefly described how to reduce the perturbed Einstein and Maxwell equations 
written in terms of the gauge-invariant variables to a set of decoupled master equa- 
tions for a single scalar variable on the 2-dimensional background spacetime for the 
vector and scalar type perturbations. For the tensor type perturbations, there is no 
electromagnetic perturbation mode and the reduction can immediately be done to 
obtain the master equation on the generic m-dimensional spacetime. In the black 
hole background case, our master equations generalize to higher dimensions in a man- 
ifestly gauge-invariant manner the well-known Regge-Wheeler-Zerilli equations for 
Schwarzschild black holes, and Moncrief's equations'^®^'*^^^ for Reissner-Nordstrom 
black holes, as well as the master equations given by Cardoso and Lemos^^) that 
include cosmological constant in 4-dimensions. We have seen that by taking Fourier 
decomposition with respect to the time coordinate, each of the master equations is 
expressed in the form of a one-dimensional self-adjoint ODE. Therefore, as guaran- 
teed by spectral theory of self-adjoint operators, the master equations can govern 
all possible perturbations that are normalizable with respect to the standard inner 
product. This is in particular important when we address the stability problem of a 
given solution, as the stability proof should be a statement concerning all physically 
acceptable (normalizable) perturbations. This is in contrast to the rotating black 
hole case, for which we have Teukolsky's equations in 4-dimensions and a similar set 
of master equations®^"^^-* for some special cases of higher dimensional Myers-Perry 
black holes, but those master equations do not reduce to a form of the self-adjoint 
eigenvalue problem, except for some special modes [c.f., Ref. 25)]. In section 7, 
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we have also briefly reviewed a similar type of master equations for gravitational 
perturbations of static black holes in generic Lovelock theory, derived recently. 

As an immediate and one of the most important applications of the master 
equations, we have examined the stability of higher dimensional static black holes 
with charge and cosmological constant in section 8. The task is to show the positivity 
of the self-adjoint operator appeared as the spatial derivative part of the master 
equations. For higher dimensional black holes, the potential term of the relevant self- 
adjoint operator is in general not positive definite, in contrast to the 4-dimensional 
case, and the stability is therefore not taken for granted in higher dimensions. To deal 
with such situations, we have developed the S'-deformation technique, in which our 
task is to find some suitable function, S, that makes the effective potential deformed 
by S positive definite. Having applied this technique we have shown that a large 
variety of static black holes in higher dimensional general relativity are stable under 
gravitational as well as electromagnetic perturbations as summarized in Table I. This 
technique has also applied to the (in)stability analysis in general Lovelock theory, as 
has just been discussed in section 8.5 above. A similar type of stability analysis using 
the S-deformation, restricted to tensor-type perturbations, has also been performed 
for static black holes in higher derivative stringy gravity.^^^ Thus, the S'-deformation 
has turned out to be a powerful tool to address stability problem. 

However, as also indicated in Table I, the stability analysis of higher dimensional 
black holes has not been completed yet even within the context of general relativity. 
For some cases, in particular, when electric charge and cosmological constant are 
involved, we have not yet been able to draw definite conclusions for the stability 
problem with respect to scalar type perturbations. This is in part because there does 
not seem to be a systematic method to find such a desirable function S that could 
apply to generic cases. In particular, for the charged black hole case, the potential 
for scalar-type perturbation admits a negative ditch in the immediate vicinity of the 
horizon, which appears to be difficult to remove by the ^-deformation. In connection 
to this, it would be interesting to note that the existence of such a negative ditch in 
the potential may have a significant infiucncc on the frequencies of the quasinormal 
modes and the graybody factor for the Hawking process, even if these black holes 
are found to be stable. Also, given a choice of S, the stability would still depend 
upon the range of the parameters characterising the black hole solution as well as 
upon the eigenvalues of the harmonic functions. Therefore our analysis using the 
iS-deformation needs to go on a case-by-case basis. A numerical analysis^^^ has found 
no indication of instability for charged AdS black holes in D = 5, 11. Therefore, 
at least for charged AdS black holes it may still be possible to analytically prove its 
stability by using the S'-deformation or other analytic methods. 

As other applications than the stability issue, the master equations can be used 
in numerical studies of black hole quasinormal modes in higher dimensions.^^^'^^^"^^^ 
It would also be interesting to consider stationary perturbations that could describe 
deformation of the event horizon, as considered in Ref. 37) for unstable charged 
de Sitter black holes. If one finds no stationary perturbation that is regular ev- 
erywhere on and outside the event horizon, then it would support the uniqueness 
property^^)'^^) of the given background black hole solution, as analysed in Paper II 
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for the vacuum black hole well as in Ref. 57) for more general cases. This may 

be interesting in particular in asymptotically AdS black hole case to find deformed 
horizon solutions [see Ref. 58) for such an analysis in 4-dinicnsions]. Also, we note 
that depending on the type of perturbations and the dimensionality, asymptotically 
AdS black holes admit a large class of boundary conditions at conformal infinity, 
other than the Dirichlet conditions considered in section 8. As has been consid- 
ered in the context of gauge/gravity correspondence and often examined within AdS 
gravity coupled to a scalar field, it would be interesting to clarify whether dif- 
ferent choice of boundary conditions leads to different consequences for the stability 
problem. 
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